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unstable transitions) are reanalyzed. A new kind of plot is presented for 
studying phase transitions in finite systems such as atomic nuclei. The im- 
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I. INTRODUCTION 

In the last few years, interest for the study of phase transitions and phase coexistence 
in atomic nuclei has been revived |l Ij in particular making use of the Interacting Boson 
Model (IBM) 0. In the chart of nuclei, three transitional regions can be distinguished 
where one observes rapid structural changes: (a) In the Nd-Sm-Gd region, one observes a 
change from spherical to well-deformed nuclei when moving from the lighter to the heavier 
isotopes; in the IBM language this is the f/(5)-577(3) transitional region; (b) in the Ru-Pd 
region, one notices that the lighter isotopes are spherical while the heavier ones indicate a 7- 
unstable character, this is the U(5)-0(6) region; (c) in the Os-Pt region, the lighter isotopes 
are well deformed while the heavier shown 7- unstable properties, this is the SU(3)-0(6) 
transitional region. Although these three transitional regions have been studied extensively 
in the framework of the IBM |§-|9| (note that in this paper we use the simplest version of the 
model, IBM-1, therefore we do not consider the transitional regions using IBM-2 |7|,[lT)|,|TT|j), 
the discussion of phase transitions has not always been treated in a proper way. In particular 
one of the weak points is how to define an appropriate control parameter in a system such 
as the atomic nucleus where this parameter is fixed (for given iV and Z) and cannot be 
controlled externally ||12|| . This problem has been recently considered in a study by Casten 



et al. §. 

Binding energy (BE) serves as an appropriate "signature" for phase transitions in nuclei 



12|jl3|| . The binding energy of the nucleus or, equivalently its mass, is one of the more 
important nuclear properties to be determined. In some cases this is the only knowledge 
about a nucleus when situated far from stability. In the last few years, the development 
of new experimental facilities has given access to very unstable nuclei and many new mass 



measurements have been obtained [|1J]. In the light of the variety of nuclear properties 
(nuclear mass, nuclear moments, . . .) that have been obtained recently, a theoretical study 
concentrating on the nuclear binding energy and, related to that, possible phase transition 
is desirable. 



Binding energies contain important information about nuclear structure. In IBM calcu- 
lations, only spectra and electromagnetic transitions used to be analyzed and the description 
of the binding energy is in most of the cases excluded. In the present article we show that it 
is convenient to explore nuclear binding energies in order to obtain a consistent description, 
especially when one is dealing with chains of isotopes. It must be noted that the study 
of binding energies can help to discriminate between apparently equivalent descriptions of 
energy spectra and electromagnetic transitions rates for a given nucleus or, even for a chain 
of nuclei. 

In section || we review the study of phase transitions in the IBM. In section [III] a new 
way of dealing with phase transitions in nuclei will be presented together with the study of 
two-neutron separation energies (S271) in regions where phase transition might happen. In 
section |TV| we study some realistic cases and, finally, in section |V|, we present our conclusions. 



II. TRANSITIONAL REGIONS AND PHASE TRANSITIONS IN THE IBM: A 
STUDY OF NUCLEAR BINDING ENERGIES 

The IBM is an algebraic model able to describe the low-lying collective states in even-even 
nuclei in terms of bosons which carry either angular momentum L = (s bosons) or angular 
momentum L = 2 (d bosons) ||. The system of bosons, for which the number of bosons 
equals half the number of valence fermions, N = n/2, interacts through a Hamiltonian 
that typically includes up to two-body interactions, is number conserving and rotationally 
invariant. The original version of the model, the IBM-1, is used in the present article. 

The IBM has been extensively used during the last decades for the study of medium 



and heavy nuclei [|15 |. One of the main challenges for the model was the description of 
transitional regions, where the structure of the nuclei can change rapidly along a chain of 
isotopes. In particular, different transitional regions have been obtained where such changes 
in the nuclear structure appear: from spherical to well deformed nuclei (Nd, Sm, and Gd 
region), from well deformed to 7-unstable nuclei (Ru and Pd region) or from spherical to 



7- unstable (Os and Pt region). In the IBM language a spherical nucleus is related to the 
U(5) limit, a well deformed nucleus is related to the 577(3) limit, while 7-unstable nuclei 
correspond to the 0(6) limit. A compact Hamiltonian for analyzing the three previous 
regions can be depicted as follows, 

H = n ( yN l_ifi d -Q-Q)+n'L-U (1) 

with < £ < 1. The coefficient of the one-body term can be rewritten as N(l — £)/£ = e^/ft, 
where e^ is the single-particle energy for the d bosons. On the other hand, fid is the d boson 
number operator and 

L = VlO(S x d) {1 \ (2) 

Q = s t d + d t s + x(rf t x d) (2) . (3) 

The symbol • represents the scalar product. In this work the scalar product of two operators 
with angular momentum L is defined as T L ■ T L = J2m(.~^-) M TlmTl-m where T LM corre- 
sponds to the M component of the operator T^. The operator 7^ m = (— l) m 7^_ m (where 7 
refers to s and d bosons) is introduced to ensure the tensorial character under spatial rota- 
tions. With this Hamiltonian one can explore the three transitional regions which correspond 
to the segments of the symmetry triangle [|l^ . 



One of the most striking facts that happen in the transitional regions is the possibility of 
observation of phase transitions. A phase transition can roughly be defined as a qualitative 
change in a given property of a system. In order to characterize the properties of the 
system one has to introduce an order parameter. Another important concept is the order 
of a transition which is the degree of the derivative for a given experimental observable 
which first experiences a discontinuity. In the atomic nucleus this experimental observable 
is the binding energy or, the related two-neutron separation energy (S2n)- Although phase 
transitions are strictly well defined in macroscopic systems and the atomic nucleus is a finite 
system, a number of studies have shown that the concept of a phase transition retains its 



validity and usefulness in small systems too [|17| , |18|,[13|,[19fl . A very useful tool in order to 



discuss phase transitions in finite system are the coherent states which, in the case of the 



IBM, are also known as intrinsic states PD,|l2,BI|. The intrinsic-state formalism provides 



a connection between the IBM and the Geometrical Collective Model ||22|| . The key point 
for establishing this connection is to consider that the dynamical behavior of the system 
can be described, to a first approximation, as arising from independent bosons moving in 
an average field [23]. The ground state of the system is a condensate, |c), of bosons that 



occupy the lowest-energy phonon state T' c , 



\c) = -4f(TjnO>, (4) 



where 



r ^ = 7TTW ( st + / 3cos 74 + ^|/ 5sin 7(4 + d- 



(5) 



Here, the variables j3 and 7, that are obtained after minimization of the ground-state energy, 
are the order parameters of the system. A spherical nucleus will have as order parameters 
(3 = and an arbitrary value for 7; a deformed nucleus is characterized by a finite value 
of (3 (in the SU(3) limit for N — ► 00 one has (3 = \/2) and 7 = (prolate nucleus) or 
7 = 7r/3 (oblate nucleus); and a 7-unstable nucleus corresponds to (3 = 1 (in the limit 
iV — > 00) and to an arbitrary value for 7. A phase transition will produce changes in the 
order parameter of the system. Using the intrinsic state we get "microscopic" information 
on the order parameter. 

First, we analyze the behavior of the binding energy and two- neutron separation energy 
in the three transitional regions. In the three cases the results are independent of the value 
of k' in eq. ([]]). For the U(5)-SU(3) transition the value of £ varies from (£7(5)) to 1 
(577(3)) with \ — iV^/2 [fTE| ; in the f/(5)-0(6) transition the value of £ varies from 
(£7(5)) to 1 (0(6)) with x = JIB]] ; finally in the transition 577(3)-0(6) the value of x varies 
from iv^/2 (577(3)) to (0(6)) with £ = 1 fllBJ. In figures 0a, |l]b, and |l|c we present the 



binding energy in the laboratory system and in the intrinsic frame (BE/ k) while in figures 
|l]a', ITJb', and IJc' we present the value of (3 in the three transitional regions for a number 



of bosons N = 10. First of all, it must be emphasized that the binding energy calculated 
in the intrinsic system is always below the exact value because it derives from a variational 
method. One can also realize that the intrinsic-state calculations show a sharp behavior 
in the phase transition regions while the diagonalization in the laboratory frame leads to 
smooth curves. In this sense the variational procedure mimics very well what happens in 
the thermodynamic limit. One of the characteristics of phase transitions in macroscopic 
systems is that they are very sharp, but in finite system the sharpness will depend on the 
number of particles. 

In figure [3]a, a first order transition shows up (discontinuity in the slope of the BE) , and 
also a discontinuity in the value of the order parameter exits. In figure |l|b, a phase transition 
also appears, but in this case it is a second-order transition (discontinuity in the curvature 
of the BE). One observes, more generally, that there exits a line of first order transitions 
for Hamiltonians with x 7^ and just one point for a second-order transition in the case 
of x = 0- The parametrization of the Hamiltonian ([I]) is very appropriate because the 
phase transition appears approximately at the same value of the control parameter (£) for 
any value of N [|J. This critical point, £ c , can be obtained easily from the expression for the 
PI 



energy given in 



as 



N 

For £ < £ c the potential has a local minimum in (3 = 0, while for £ > £ c this local minimum 
becomes a local maximum. It is clear from this expression that for a large number of 
bosons, £ c m 0.2. Finally, in figure |TJc no phase transition appears and both calculations, 
exact diagonalization and intrinsic-state, lead to a smooth transitions. 

In a next step, we analyze the behavior of the two- neutron separation energy. The 
definition that will be used along this paper is, 

S 2n (N) = BE(N)-BE(N-l). (7) 

In the calculation of S2« two different nuclei are involved, which makes necessary to impose 
an assumption in order to relate the Hamiltonians for the nuclei with N and N — 1 bosons. 
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Here, we consider the simplest (not unrealistic) case which implies the same Hamiltonian in 
both nuclei. In figure |2| we study the same transitions as in figure [I] and essentially the same 
remarks hold. Note in this case the crossing of intrinsic-state and laboratory calculations, 
which is allowed. 

So, in this section we have reviewed the study of phase transitions in atomic nuclei, in 
particular using the IBM, and have shown it is straightforward and clear from a theoretical 
point of view. In this analysis, however, there is a weak point: the control parameter £ 
(x) is not a genuine one, because it cannot be modified externally. The control parameter 
only changes when one is moving along a chain of nuclei and also assumes a change in the 
number of bosons. These two problems have already been remarked and alternative solutions 
have been proposed |J. In the next sections, we propose a new strategy for attacking this 
problem. 

III. HOW TO IDENTIFY PHASE TRANSITIONS 

A. The new diagrams "binding energy-number of bosons-control parameter" 

A consistent treatment of phase transitions in nuclei must go beyond the study carried 
out in section [II| because one has to treat simultaneously the control parameter £ (%) and 
the number of bosons N. It must be emphasized that N is not a control parameter because 
it is a discrete variable and it is fixed in each individual nucleus. Strictly speaking £ (x) is 
neither a genuine control parameter, because its value is also fixed in each given nucleus, 
nor can it be controlled externally. Only the appropriate interplay of N with £ (x) can lead 
to a proper definition of an order parameter, although with discrete values. 

The appropriate way for treating phase transitions and transitional regions is to plot in 
the same figure the value of the binding energy versus the control parameter and this for 
different values of N. In such a plot, a transition will develop through changes in the control 
parameter, but at the same time going through curves with a different number of bosons. 



In order to illustrate this new procedure, we discuss the transitional regions U(5)-SU(3), 
U(5)-0(6) and, SU(3)-0(6) in figures ||a, |3|b, and ||c respectively. In these figures, only the 
laboratory results (IBM diagonalization) are presented. As an example, trajectories for real 
nuclei are also plotted in each figure. A procedure in order to obtain the value of the control 
parameter, for each nucleus, will be explained in section |TV[ The more interesting cases 
appear in figures ||a and ||b because there, phase transitions indeed happen. So, looking at 
these figures, one can easily see if a given nucleus is situated in the spherical region (£ near 
to 0), in the deformed region (£ around 1) or at the critical point. 

One might be tempted to consider similar plots for the two-neutron separation energies. 
In this case one must be very careful because the Hamiltonian for N and N — 1 bosons 
are different. That means that S2 n is not a function of a single control parameter, £ (%), 
anymore, but becomes a function of two control parameters, one for N and another one for 
N-l. 

The essential difference between phase transitions in nuclei compared to macroscopic 
systems is the existence of a control parameter in the macroscopic case, such as temperature 
or pressure, that can be controlled externally. In the case of atomic nuclei, the only way 
for changing the value of the control parameter is to move along a chain of nuclei. It is 
clear that the variation of the control parameter is thereby fixed and cannot be modified 
externally. 

B. Crossing the phase transition region: anomalies in the two-neutron separation 

energies 

In previous sections we have studied the behavior of the binding energy as a function 
of the control parameter, i.e. as a function of the Hamiltonian, in the three transitional 
regions of the IBM. Due to the specific characteristics of phase transitions in atomic nuclei, 
it is not a priori clear how different observables will behave when crossing a critical value 
through a chain of isotopes. More particular, one can observe the two-neutron separation 
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energies as a function of the nucleus in the chain that one considers. S 2n is indeed a very- 
appropriate observable to be analyzed because it contains important information about the 
nuclear structure, in particular on the ground state. 

Before starting the present analysis we emphasize the fact that there exits an extra 
contribution to the binding energy not yet included when using the Hamiltonian of equation 
(I), ie. 

E' = E + AN + ^N(N - 1) (8) 

Those extra terms, deriving from the linear and quadratic U(6) Casimir invariant, take into 
account the bulk contribution of the nuclear interaction. Including those extra terms in S 2n , 
one obtains a linear contribution to be added to the IBM calculation, resulting 

S 2n (N) = (A- B/2) + BN + BE IB m(N) - BE IBM (N - 1). (9) 

When the linear dependence is excluded from S 2ra , we will refer to that result as S' 2n . It 
must be noted that the specific nuclear structure IBM contribution in the three dynamical 
limits (U(5), SU(3), and 0(6)) is also linear in the boson number P,|25| and in this context 
only a particular change in the internal structure of the nuclei, as happens during the phase 
transitions, can perturb this linear dependence. 

In order to be sure one encounters a phase transition or passes through a specific tran- 
sitional region, one has to establish the value of the control parameter for each nucleus. In 
the following we will concentrate on chains of isotopes (fixed number of protons), being the 
nuclei characterized by a given number of bosons iV or by the number of neutron pairs N„. 
As a consequence one has to determine the functional relation £ = £(iV) (x = x(N)) or, 
in terms of the number of neutron pairs, the relation £ = £(N V ) (x = x{N u ))- Although 
the latter functional relation must be determined from experimental data, one can use a 
simple function and see if it is possible to obtain some physical insight on the structure of 
S 2n . For the present study we fix a chain of isotopes with 5 protons pairs, N n = 5, and 
a variable number of neutrons pairs, ranging from N u = to N v = 10. Two functional 



dependences will be used, firstly linear and secondly quadratic. Next, we show the different 
parametrizations for the different transition regions, 

• U(5)-SU(3) and C/(5)-0(6) transitional regions. 

£ = 0.099iV„ + 0.01, (10) 

£ = 0.0099A^ + 0.01, (11) 



SU(3)-0(6) transitional region. 



X = -%N. (12) 



Using parametrizations (|10|,|TT]), £ ranges from 0.01 to 1. On the other hand, using 



parametrization (|12| , |T3|) , x ranges from to — vf/2. In order to establish a realistic 
parametrization, we start from the empirical observation that in the U(5)-SU(3) transi- 
tion, the system passes from the U(5i) to the SU(3) limit when the number of bosons is 
increasing |]15|| . Similarly, the transition U(5)-0(6) implies that the system goes from U(5) 
to 0(6) when the number of bosons is increasing. Finally, the transition SU(3)-0(6) implies 
the system passes from 577(3) to 0(6) when the number of bosons is decreasing. 

In figures |]a and||b, S2 n and S^n are shown, respectively, for the U(5)-SU(3) transitional 
region using for £ the parametrizations (|lOl , |TTl) . The Hamiltonian (HD has been used with 
X = — v^7/2 fl6l1 . For convenience, a linear dependence equal to S 1 ^ / k = 200 — 20^ has 
been chosen as a reference value in order to obtain a realistic behavior in S2 n - In these 
figures, an anomaly in the linear dependence of S 2n and S' 2n , right at the place where 
the phase transition happens, can be appreciated. The phase transition point corresponds 
approximately to £ ~ 0.2 (see figure |l]a and 0a') and taking into account equations (0) and 



11]) results a value for N u , N u ps 2 and N v ~ 4 for the linear and quadratic dependence, 



respectively. We point out that in case a quadratic variation in N u for the control parameter 
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is selected, the anomaly in S 2n is more pronounced than for a linear variation. In figures [|a 
and |b, S 2n and S' 2n are shown, respectively, for the U(5)-0(6) transitional region using for 
£ the same parametrization as in the previous case. The Hamiltonian (|IJ) has been used with 
X = [[LJJ]. In this case a linear dependence equal to S 1 ^ / 'k = 200 — lOiVj, has been used 
to determine the S 2n values. Again, a non-linear behavior appears at the point of the phase 
transition. The phase transition point corresponds approximately to £ ~ 0.2 (see figure |l]b 
and |l|b') and taking into account equations (|12|) and flilf ) results a value for N v , N v « 2 and 
N v m 4 for the linear and quadratic dependence, respectively. In this case, a linear variation 
in the control parameter leads to an almost linear dependence in S2 n and S' 2n . Only a control 
parameter which depends quadratically on N v produces a kink in the two-neutron separation 
energy. Finally, in figures ||a and |]b, S 2n and S' 2n are represented, respectively, for the SU(3)- 
0(6) transitional region using for \ the parametrizations (0,0). The Hamiltonian ([!]) has 
been used with £ = 1 [ 16| . In this case the following linear dependence S 1 ^ / K = 200 — 20^ 
has been chosen as a reference when calculating S2 n - In this transitional region no phase 
transition is observed and a smooth curvature result except when one is approaching the 
SU(3) limit and is using a quadratic dependence in N v for \ (0)- 

Along this section we have learned to treat in a proper way binding energies (two-neutron 
separation energies) in transitional regions where phase transitions might appear. We have 
seen that in these regions one can get deviations from the linear tendency of S 2n - Next step 
will be to find experimental examples where those anomalies appear and explain them using 
the IBM. 

IV. REALISTIC CALCULATIONS 

Within the framework of the IBM, the energy spectra and transition rates of many 
medium-mass and heavy nuclei have been analyzed very successfully. However, in most of 
these studies [|TJ|] the binding energies have been ignored. One of the reasons for such a 
neglect stems (partly) from the fact that the early IBM studies (carried out in the 80 's) 
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|7| JT0|JTT[| gave a rather good reproduction of the general trends of the nuclear binding energy 



in different mass regions. So, it was tacitly assumed that realistic IBM calculations would 
also reproduce (in a natural way) the binding energy values. In the mean time, with a 
large increase in the precision for the mass measurement p?| , p8| , p^| , much improved data 



have been established [ESJ. Moreover, it must be noted that, in order to carry out a proper 



comparison between calculated and experimental binding energies (or S2 n data), one cannot 
concentrate on single or a few nuclei. Large chain of isotopes have to be considered and this 
complicates considerably performing consistent IBM calculations. 

In this section we will determine the experimental relations between the control param- 
eter £ (x) and the number of neutron pairs N v through explicit fits of energy spectra for 
selected nuclei in the three transitional regions. Starting from these relations we are able 



to represent a specific chain of isotopes on the diagrams shown in section |T H A| . We subse- 
quently show that the experimental binding energies (or two-neutron separation energies) 
are not always well reproduced although the spectra do. We analyze chains of isotopes in 
the transitional regions (i) first using the more standard set of parameters for the Hamilto- 
nian H and (ii), secondly, we compare with an alternative parametrization (see discussion 
below). We do not aim a perfect fit of the experimental data, but we intend to emphasize 
the influence of choosing various parametrizations on the binding energies. 
Next we study the three transitional regions: 

• The U(5)-SU(3) transitional region 

There exit nuclei that exhibit a drastic change in the internal structure when one 
is passing their isotopes series, going from a vibrational into a rotational behavior. 
Typical examples are situated in the rare-earth region, such as the Nd, Sm, and Gd 
nuclei. The Sm nuclei have been extensively studied in Ref. |J and a very good 
agreement with experimental data for S2n has been obtained. In figure [3]a, we plot 
the position of the different Sm isotopes and one can clearly see how they cross the 
critical point, in particular 152 Sm is situated almost exactly where the phase transition 
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happens. The cases of Nd and Gd are very similar. Here, we will only focus to the 
case of Gd. 

In this region we assume a value x — ~ V*/2 because some Gd isotopes clearly exhibit 
the character of the 577(3) dynamical symmetry J5|. This assumption was very suc- 
cessful in describing the Sm nuclei, which form neighboring nuclei. With this ansatzwe 
have carried out a fit to the energy spectra of the Gd isotopes (from 150 Gd until 162 Gd) 
using the Hamiltonian ([I]). In order to obtain the parameters, we tried to reproduce 
the rotational structure of the ground- and 7-band as well as possible. The 150 Gd 
nucleus still shows a vibrational structure while 156_162 Gd are considered as rather 
good SU(3) examples. The parameters for 150 Gd to 162 Gd, respectively, are shown in 
table H (upper part). In figure 0a we compare the theoretical and the experimental 
energies for the low-lying levels using this set of parameters. Those values are used 
to position the Gd isotopes in the phase transition diagram of figure ^a, where 150 Gd 
corresponds to iV = 9 while 156 Gd corresponds to N = 12. Note that 158_162 Gd are 
not plotted; they also corresponds to £ = 1. From figure [3|a one can observe a sudden 
transition in the Gd isotopes from a vibrational regime into the 577(3) limit. Up to 
now, we have obtained a good description of the Gd energy spectra but we have to 
consider the experimental two-neutron separation energies too. This comparison is 
carried out in figure [| (the present calculation corresponds to Theo.-a), using a linear 
contribution S%™ = 15 — 0.252N U (MeV) as a background. One observes a very bad 
agreement which points to the fact that the use of x — ~ V7/2 is not appropriate. A 



different approach to fix the Hamiltonian is proposed in [BB] . There, it is assumed that 
in well deformed nuclei one can keep the value of k constant around k = 30 keV and 
smoothly change the values of x an d £d along a series of isotopes. In the present case, 
the appropriate values for \ and k are % = — 0.6 and k = 19.2 keV. On the other hand 
the values of £ (using N(l — £)/£ = e d /n, see equation ([I])) are shown in table | (lower 
part) in going on 146 Gd until 162 Gd (note that in this case, a few extra isotopes are 
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also considered). In figure [7|b we again compare the theoretical and the experimental 
energies for the low-lying levels corresponding to the new Hamiltonian. Again, we 
compare the calculated S2n values with the data and, as can be seen in figure [8] (this 
calculation corresponds to Theo.-b), using a linear contribution S 1 ^ = 16.75 — 0.607^ 
(MeV) as a background, a very good agreement is obtained. Note that the present 
results for £ = £(iV) cannot be represented in figure ||a due to the new value of \i 
since figure |||a corresponds to x = ~ v7/2. It is worth noting the strong similarities 
between figure || and figure f§b. 

The U(5)-0(6) transitional region 

Clear-cut examples of nuclei in which the structure changes from spherical to 7- 
unstable are obtained in the Ru and Pd region. Both series of isotopes have already 
been studied in ||, showing that the lighter isotopes exhibit a vibrational structure 
while the heavier ones present a 7-unstable behavior. In reference || the general trends 
have been analyzed but it was shown that difficulties appear when trying to obtain 
a consistent description for several observables at the same time. Due to the strong 
similarities between both nuclei, only the study of Pd nuclei will be presented here. 

A first assumption will be to consider x — 0, which is suggested in reference || too. 
Using the Hamiltonian ([[]) we fit the energy spectra of 100_112 Pd. Nc> other isotopes 
are fitted because only for those the energy ratio E(4f)/E(2i) changes between 2 
(U(5)) and 2.5 (0(6)). Note that this energy ratio will be the main observable to be 
reproduced during the fits. We have taken a fixed single particle energy at e^ = 65 keV 
and the parameters that we have obtained are shown in table || (upper part). Note 
that no rotational term is necessary. In figure |9]a we compare the theoretical and 
the experimental low-lying levels for this set of parameters. Those values are used 
to position the Pd isotopes on figure [3]b, where 100 Pd corresponds to iV = 4 and 
112 Pd corresponds to N = 10. One can clearly see that the Pd isotopes cover only 
half the range between the U{5) and the 0(6) limit. In figure one obtains a poor 
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description when comparing S2n (this calculation corresponds to Theo.-a), with a linear 
contribution S 1 ^ = 22.27 — 1.295./V,, (MeV) as background, while the energy spectra 
are well reproduced. Again, this suggests a change in the value of x- I n the following, 
we will use \ = —0.3, a value that produces a realistic description of the energy 
spectra. In this case one can fix the single particle energy to e^ = 70 keV. The 
obtained values of (k,£) from 100 Pd until 114 Pd are shown in table fl| (lower part). In 
figure ||b we compare again the theoretical and the experimental low-lying levels for 
the new Hamiltonian. Now, however, one obtains good agreement when comparing 
experimental and theoretical S271 values (this calculation corresponds to Theo.-b), with 
a linear contribution Si^ = 20.46 — 0.835iV„ (MeV) as a background, as can be seen 
in figure [H]. The new results cannot be represented on figure |3|b because x^O- 

The S77(3)-0(6) transitional region 

The most clear-cut examples of nuclei situated in this transitional region are the Os 
and Pt nuclei. The lighter isotopes exhibit rotational structures while the heavier ones 
come very close to the 7- unstable limit. This region has already been studied in || in 
the framework of the IBM-1 but using a different parametrization for the Hamiltonian. 
Due to the similarities among both series of isotopes, only the analysis of the Pt nuclei 
will be presented here. 

In order to make a global fit we take into account the fact that 196 Pt is one of the 
best examples of a 7-unstable nucleus and, as a consequence, x = and £ = 0. 
The obtained parameters («, x, k?) for 184 Pt until 196 Pt are shown in table |t| (upper 
part). In figure |ll|a we compare the theoretical and the experimental low-lying levels 
for these parameters. If one places these values in figure [3]c, where 196 Pt corresponds 
to N = 6 while 184 Pt corresponds to N = 12, it is clear that the Pt isotopes always 
stay very close to the 0(6) limit. Finally, in figure O, we compare the theoretical 
and experimental S2 n values (this calculation corresponds to Theo.-a), with a linear 
contribution S^n — 12.66 — 0.545A^, (MeV) as a background. Again, we obtain a poor 
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agreement in the description of the experimental values. In order to treat this mass 



region with a slightly different Hamiltonian, one can follow again Ref. ||26| . There, it 
is suggested to include a single-particle energy term. In the case of the Pt nuclei the 
more appropriate value is e^ = 25.7 keV. The value of k is also fixed to K — 33.5 keV 
and k' = 15.2 keV. The different values of the parameters in the Hamiltonian for 184 Pt 
until 196 Pt are shown in table |l| (lower part). In figure |TTl b we know compare the 



theoretical and the experimental low-lying levels for the new Hamiltonian. Using this 
new parametrization, when comparing the S2 n values (this calculation corresponds to 
Theo.-b), with a linear contribution S 1 ^ = 11.85 — 0.649^ (MeV) as a background, a 



good agreement is obtained as can be seen in figure [12]. These new results cannot be 
represented in figure |]c because now e^ ^ 0. 

In this section we have presented two different parameter sets obtained from fitting the 
energies of the low- lying states of different isotopes of Gd, Pd, and Pt. A very striking fact 
is, on one side, the good description of the energy spectra for both sets of parameters, and 
on the other side, the very different results obtained for the two-neutron separation energy 
using these different parametrizations. The results obtained here were totally unexpected 
because it was assumed tacitly that a good description of energy spectra would imply a good 
description for the binding energies too. In order to illustrate this problem, we consider the 
case of 148 Sm which has been studied by different authors P,p6|. Although in all the cases 
the description of the energy spectrum is satisfactory, when calculating the binding energy 
for the proposed Hamiltonians, one obtains the results 0.939 MeV (Ref ||), 1.571 MeV 
(Ref f26[) and 2.253 MeV (also in Ref p6|). All this suggests that a consistent study of 
atomic nuclei, within the IBM, requires the inclusion of an analysis of the binding energy. 

V. SUMMARY AND CONCLUSIONS 

In the present article, we have analyzed transitional regions and phase transitions in 
the framework of the IBM. A new kind of diagram (see figure |]) for dealing with phase 
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transitions has been presented. In this diagram we plot the binding energy of the system 
versus the control parameter (Hamiltonian) for a range in the number of bosons. This 
procedure turns out to be a particularly convenient method when working in finite systems 
without a true control parameter. 

The main conclusions are twofold: 
(i) We have shown that, although the crossing through the phase transition points produces a 
deviation in the S2 ra values from the overall linear background dependence, the way in which, 
in the given chain of isotopes, this critical point is crossed proves extremely important in 
deriving clear-cut results. If e.g. the relation between a control parameter and the isotope 
location {e.g. boson number) is linear, one does not observe dramatic changes in S2 n in any of 
the three transitional regions (U(5)-SU(3), C/(5)-0(6), and S77(3)-0(6)). On the contrary, 
for a quadratic relation, clear changes can appear in the S2 n behavior at the crossing point, 
(ii) We have observed that in a consistent study of long chains of isotopes, one has to treat 
the ground-state (through its binding energy) on equal level with the excited states (relative 
energy spectrum). It turns out that parameters producing very similar energy spectra can 
still result in important differences in the ground-state binding energy (order of 1 MeV). 
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FIG. 1. Binding energies and (5 as a function of the control parameter, in the transitional 
regions (a) U(5)-SU(3), (b) U(5)-0(6), and (c) SU(3)-0(6). Dashed lines correspond to laboratory 
calculations while full lines correspond to intrinsic state results. Number of bosons N = 10. 
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FIG. 2. Two-neutron separation energies (S2n) a s a function of the control parameter, in the 
transitional regions (a) U(5)-SU(3), (b) C/(5)-0(6), and (c) SU(3)-0(Q). Dashed lines correspond 
to laboratory calculations while full lines correspond to intrinsic state results. Number of bosons 

N = 10. 
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FIG. 3. Binding energies in the laboratory system for different numbers of bosons, N, as a 
function of the control parameter, in the transitional regions (a) U(5)-SU(3), (b) U(5)-0(6), and 
(c) SU(3)-0(6). Full circles and open squares correspond to the theoretical positions of different 
chains of isotopes. 
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FIG. 4. Schematic plot of (a) S2n and (b) S\ n in the U(5)-SU(3) transitional region. Full 
lines corresponds to a linear variation of £ with respect to iVj, while dashed lines corresponds to a 
quadratic dependence. 
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FIG. 5. Schematic plot of (a) S2n and (b) S\ n in the U(5)-0(6) transitional region. Full lines 
corresponds to a linear variation of £ with respect to N v while dashed lines corresponds to a 
quadratic dependence. 
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FIG. 6. Schematic plot of (a) S2n and (b) S\ n in the SU(3)-0(6) transitional region. Full 
lines corresponds to a linear variation of x with respect to N v while dashed lines corresponds to a 
quadratic dependence. 
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FIG. 7. Comparison between theoretical and experimental energy levels for Gd isotopes. Full 
lines correspond to IBM calculations, while symbols connected by dashed lines correspond to 
experimental data. In panels (a) and (b) different sets of parameters for the IBM calculations are 
used, see text. 
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FIG. 9. Same caption as in figure [7| but for Pd isotopes. 



28 



> 



CO 



20.5 

20 
19.5 

19 
18.5 

18 
17.5 

17 
16.5 

16 
15.5 

15 
14.5 

14 
13.5 

13 
12.5 







I 




- 


\ * 






- 


♦-- 


w l_A|J. 

-♦ Theo.-a 


- 


^^ '■■•• 
v '"•■ 
\\ \ 




■- - 


--■ Theo.-b 


- 








- 


- 






; 


- 




*"■■ «x 


i , 


- 




'■'-♦"' 


\^j 






i 





5 
N. 



FIG. 10. Same caption as in figure but for Pd isotopes. 
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FIG. 11. Same caption as in figure ^ but for Pt isotopes. 
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FIG. 12. Same caption as in figure but for Pt isotopes. 
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TABLES 
TABLE I. Parameters of the Hamiltonian for Gd isotopes (N w = 7). Upper part of the table 

is referred as Theo.-a and lower part as Theo.-b (as used in the text and in the figures). 

Theo.-a 



A 


150 


152 


154 


156 


158 


160 


162 


N u 


2 


3 


4 


5 


6 


7 


8 


K 


15.4 


15.4 


15.4 


15.4 


14.8 


11.3 


9.1 


£ 


0.139 


0.192 


0.287 


1 


1 


1 


1 


k' 


9.0 


9.0 


9.0 


9.0 


7.7 


8.3 


8.6 



K, k' in keV and £ dimensionless, x = ~~ n/7/2. 
Theo-b 



A 


146 



148 
1 


150 

2 


152 
3 


154 

4 


156 

5 


158 
6 


160 

7 


162 

8 


£ 


0.60 


0.137 


0.166 


0.236 


0.373 


0.535 


0.625 


0.658 


0.724 



k = 19.2 keV, k 1 — and £ dimensionless, x — —0.6. 
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TABLE II. Parameters of the Hamiltonian for Pd isotopes (N^ = 2). Upper part of the table 
is referred as Theo.-a and lower part as Theo.-b (as used in the text and in the figures). 
Theo.-a 



A 


100 


102 


104 


106 


108 


110 


112 


N u 


2 


3 


4 


5 


6 


7 


8 


K 


20.0 


42.0 


52.0 


49.0 


47.0 


52.0 


70.0 


£ 


0.110 


0.244 


0.324 


0.345 


0.366 


0.419 


0.52 



k in keV, k' = and £ dimensionless, x = 0- 
Theo-b 



A 


100 


102 


104 


106 


108 


110 


112 


114 


N u 


2 


3 


4 


5 


6 


7 


8 


9 


K 


22.0 


44.0 


50.0 


44.0 


40.0 


37.0 


37.0 


33.0 


C 


0.112 


0.239 


0.300 


0.306 


0.314 


0.322 


0.346 


0.342 



k in keV, k' — and £ dimensionless, x = —0.3. 
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TABLE III. Parameters of the Hamiltonian for Pt isotopes (N n = 2). Upper part of the table 
is referred as Theo.-a and lower part as Theo.-b (as used in the text and in the figures). 
Theo.-a 



A 


184 


186 


188 


190 


192 


194 


196 


N u 


10 


9 


8 


7 


6 


5 


4 


K 


43.0 


44.0 


44.0 


47.0 


60.0 


60.0 


63.0 


X 


-0.115 


-0.110 


-0.080 


-0.055 


-0.049 


-0.050 





k' 


4.2 


0.8 


17.6 


19.0 


11.0 


11.0 


11.0 



k and k' in keV and £ dimensionless, Ed = 0. 
Theo-b 



A 


184 


186 


188 


190 


192 


194 


196 


N u 


10 


9 


8 


7 


6 


5 


4 


K 


33.5 


33.5 


33.5 


33.5 


33.5 


33.5 


33.5 


X 


-0.25 


-0.20 


-0.10 














k' 


15.2 


15.2 


15.2 


15.2 


15.2 


15.2 


15.2 



k and k' in keV and £ dimensionless, E& = 25.7 keV. 
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